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Abstract 

We analyze the 1496 days of SuperKamiokande data to put limits on the v e , i> e , 
+ v T and + v T components of the diffuse flux of supernova neutrinos, in 
different energy intervals and for different neutrino energy spectra. By considering 
the presence of only one component at a time, we find the following bounds at 
90% C.L. and for neutrino energy E > 19.3 MeV: < 73.3 - 154 cm _2 s _1 , 
$ Pe < 1.4 - 1.9 cm-V 1 , < (1.0 - 1.4) • 10 3 cm-V 1 and $ Pm+Pt < (1.3 - 

1.8) • 10 3 cm~ 2 s _1 , where the intervals account for varying the neutrino spectrum. In 
the interval E = 22.9 — 36.9 MeV, we find $> Ue < 39 — 54 cm _2 s _1 , which improves 
on the existing limit from SNO in the same energy window. Our results for + v T 
and + v T improve by about four orders of magnitude over the previous best 
constraints from LSD. 
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1 Introduction 

The search for the diffuse flux of neutrinos from all supernovae represents a 
new frontier of neutrino astrophysics. Just like ultra-high energy neutrinos, 
the diffuse supernova neutrinos originate, for a good part, at cosmological 
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distances. Therefore, they could give unique insights into the history of the 
universe, and specifically into the cosmological evolutions of the supernova 
rate and of the star formation rate, the physics of the first stars, etc.. This 
possibility, to do cosmology with neutrino data, is new. In addition to this, 
the diffuse flux will be a new test of the physics of core collapse supernovae 
and of neutrino propagation inside them, which will be especially important 
if a galactic supernova, a very rare event, does not occur in the next ten years 
or so. 

Until now, the searches for the diffuse supernova neutrino flux (DSNz/F) have 
turned out negative, and upper limits were put. The strongest limit comes 
from the dominant detection mode - inverse beta decay, v e +p — > n + e + 
- in the largest water Cerenkov detector available, SuperKamiokande (SK), 
above the threshold of 19.3 MeV of neutrino energy [1]. The limit on the z/ e 
component of the flux reads: 

$ Pe (£ > 19.3 MeV) < 1.2 cnrV 1 at 90%C.L. , (1) 

and is valid for certain energy spectra of the neutrino flux. This limit ap- 
proaches the theoretical predictions [2,3,4,5,6,7,8,9,10,11,12,13,14,15,16], thus 
arising the expectation that the DSNz/F may be seen in the near future, and 
triggering several theoretical studies on the subject (see e.g. the review [17] 
and references therein). 

Bounds on other neutrino species in the DSNz/F are weaker than the v e limit 
(1), due to the fact that inverse beta decay largely dominates over other de- 
tection processes in water, and that non-water detectors have smaller volumes 
than SK. 

Some attention has been devoted to the v e component of the DSNz/F. The 
strongest bound on this is indirect [18]: 

$ Ve (E > 19.3 MeV) < 5.5 cnrV 1 at ~ 98%C.L. . (2) 

It is found by converting the SK limit on z/ e , eq. (1), into a result for z/ e , using 
neutrino oscillations and the strong similarity between the fluxes of muon and 
tau neutrinos and antineutrinos produced inside a supernova. 

Among the direct limits on z/ e , the strongest is from the search for charged 
current (CC) scattering on deuterium (z/ e +d — > p + p + e~) at the Sudbury 
Neutrino Observatory (SNO), in the interval 22.9 — 36.9 MeV of neutrino 
energy [19]. The result is: 

$„ e (22.9 < E/MeV < 36.9) < 61 - 93 cm" 2 s _1 at 90%C.L. , (3) 

depending on the neutrino energy spectrum. 
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The old bounds from LSD [20] are still the best for the non-electron compo- 
nents of the flux, + v T [y x from now on) and + v T (y x from now on). 
They are obtained from searches of neutrino scattering on 12 C, and read: 



In this paper we examine the SK data published in [1] and use them to put 
limits on all the distinct components of the DSNz/F: u e , 9 e , v x and v x . 

For 9 e , we elaborate on the result (1), by generalizing the SK analysis to a 
wider range of neutrino spectra, motivated by the recent progress in the theory. 
Our result is that the bound (1) could change by up to ~50% depending on 
the neutrino spectrum. 

For the z/ e , v x and v x components of the flux, the idea is to look for events 
due to CC scattering of v e on oxygen, or due to elastic scattering on elec- 
trons of neutrinos and antineutrinos of any flavor, on top of the irremovable 
background. These events are indistinguishable from inverse beta decay, since 
they have the same signature in the detector, the Cerenkov light produced by 
the outgoing lepton. One expects that, with respect to the SNO search for v e , 
the larger SK volume could overcompensate the disadvantage of the smaller 
detection cross section, resulting in a competitive limit. A similar argument 
makes us expect a strong improvement on the LSD limits for the non-electron 
species, Eq. (4). Our numerical analysis confirms these intuitions. 

The paper is structured as follows. In sec. 2 we give generalities on the DSNz/F 
and discuss its energy spectrum. Sec. 3 is a review of the SK data and of our 
analysis of them. In sec. 4 we present results. Summary and discussion follow 
in sec. 5. 



2 Generalities 

2.1 Neutrinos from core collapse and the diffuse flux 

A core collapse supernova is an extremely powerful neutrino source, releasing 
about 3-10 53 ergs of energy within ~10 seconds in neutrinos and antineutrinos 
of all flavors in similar amounts. The energy spectrum of these neutrinos at 
the production point in the star is roughly thermal, and can be described by 
a power law times an exponential [21]. The neutrinos have average energies 
in the range of 10 - 20 MeV, with the muon and tau species having harder 



$„ B (20 < E/MeV < 100) < 3 • 10 7 cnrV 1 
®p x (20 < E/MeV < 100) < 3.3 • 10 7 cnrV 



i 



at 90%C.L. , 
at 90%C.L. . 



(4) 
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spectrum than the electron ones due to their weaker (neutral current only) 
coupling to matter. 



Inside the star, the neutrinos and antineutrinos undergo either partial or to- 
tal flavor conversion depending on the mixing angle #13 and on the mass 
hierarchy of the neutrino mass spectrum. The hierarchy is defined as nor- 
mal (inverted) if the third mass eigenstate - the one whose electron com- 
ponent is sin 2 9 13 - is the heaviest (lightest). The net conversion probability 
is due to a combination of effects of neutrino-neutrino coherent scattering 
[22,23,24,25,26,27,28,29,30,31,32,33,34,35,36] and of matter-driven resonant 
conversion (see e.g., [37,38,39]). The neutrino-neutrino effects swap the en- 
ergy spectra of the electron and non-electron components of the flux for the 
inverted mass hierarchy and above a certain critical energy, E c . Typical critical 
energies are below ~10 MeV (see e.g. [29]), with lower values for antineutrinos 
compared to neutrinos. The matter-driven conversion occurs more externally 
in the star and induces a further permutation of the energy spectra of the 
fluxes in the different flavors that enter the resonances. 

The diffuse flux of supernova neutrinos in a detector is simply the sum of 
the contributions from the individual stars. This sum can be expressed as 
an integral involving the cosmological rate of supernovae, Rsn(z), which is 
a function of the redshift z and is defined as the number of supernovae in 
the unit of comoving volume in the unit time. Starting with the present value, 
Rsn(0) ~ 0(1CT 4 ) Mpc~ 3 yr _1 , the rate increases with z and flattens at z > 1 
(see e.g. [40]). In terms of Rsn, the DSNz/F in a detector at Earth, differential 
in energy, surface and time, is given by: 

c /■*»«« . x dN det (E') dz 
HqJo dE> Jo (1 + z )3 + nA 



_ tt m (l + z) 3 + tt A 

(see e.g. [13]), where dN det (E')/dE' is the contribution of an individual su- 
pernova, inclusive of neutrino oscillations and of the redshift of energy, E 1 = 
E(l + z). Q m and Q\ are the fractions of the cosmic energy density in matter 
and dark energy respectively; c is the speed of light and H Q is the Hubble 
constant. 

Examples of calculated energy spectra for the V e component of the DSNz/F 
are shown in fig. 1. We see that the flux peaks at E rnax = 3 — 7 MeV, where it 
reaches typical values of 0(1) cm" 2 s _1 MeV _1 , and falls rapidly at E 3> E max . 
It has been shown [43] that in this high energy regime, which is relevant for 
our analysis, the flux is approximated very well (~5% accuracy or better) by 
an exponential: 

<j)(E) ~0oe- E/£o , (6) 

where the energy E is close to the energy at the peak of the spectrum, E ~ 
Emax- The total flux above a realistic detection threshold of 10-20 MeV is 
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Fig. 1. Examples of energy spectra of the v e component of the DSNi/F from the liter- 
ature, specifically from the Garching (KRJ), Lawrence Livermore (LL) and Arizona 
(TBP) models, with oscillations (see also Table 1) [21,41,42,13]. The background 
due to atmospheric v e is shown as well. This figure is taken from ref. [13], with 
permission. 

larger for larger Eq, and is proportional to the normalization of the supernova 
rate, Rsn- The large uncertainties on these two quantities translates into a 
large uncertainty on the DSNz/F. 

In Table 1 we illustrate this by summarizing several published predictions of 
the v e component of the DSNz/F. The Table gives the values of Eq found by 
fitting the energy spectra with the simple exponential (6), and the predicted v e 
flux above the SK threshold. In both columns we see large differences, due to 
the different input quantities and constraints that are used in calculating the 
DSNz/F. For example, to describe the DSNz/F energy spectrum some authors 
use neutrino spectra that result from calculations of neutrino transport (e.g., 
[21]), while others use spectra that fit the SN1987A data [15] or simply adopt 
thermal spectra with temperatures varying within indicative intervals [44]. 
This diversity of approaches results in the variation Eq = 3.8 — 8.5 MeV, which 
we use in our analysis^ An even larger variation is seen in the integrated flux, 

1 Even though the interval of Eq given here refers to electron antineutrinos, we 
use it for other neutrino species as well, as a representative range, and for easier 
comparison between different channels. 
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Table 1 



Summary of existing calculations of the i> e component of the DSNzaF. For each, we 
give the value of Eq for which a simple exponential spectrum, Eq. (6), best fits 
the neutrino spectrum in the high energy regime (E <; 10 MeV). We also give the 
integrated flux above the SK energy threshold. The models in the first six rows are 
those considered in the SK analysis [1]. In several cases, the values quoted in the 
table have been inferred from graphics in the original references. The fluxes in ref. 
[13] are calculated up to a normalization factor, which is estimated to be of order 
unity. 

which ranges from ~ 0.05 cm~ 2 s _1 to values exceeding the SK limit, eq. (1). 

In closing this section, we would like to point out that the effects of neutrino- 
neutrino scattering on flavor conversion have been studied in detail only re- 
cently, and therefore they were not included in the literature we have refer- 
enced for the spectral shape, eq. (6), and for the results in figure 1 and in 
Table 1. While the inclusion of these effects is certainly in the agenda for the 
near future (see the initial study in [36]), we believe that the literature we 
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have quoted is still a good description for the purpose of this paper. Indeed, 
older results are still valid for the normal mass hierarchy, where neutrino- 
neutrino effects are negligible. Moreover, they can reproduce (with effective 
parameters) the case of inverted mass hierarchy with neutrino-neutrino scat- 
tering above the swap energy E c . This is sufficient for our analysis, since E c is 
typically below the energy windows of interest here, E ^ 10 MeV (see Table 
6; a possible exception is LENA at Wellington, as the table shows). 



3 The data and the analysis 



3. 1 The data and their interpretation 



We have used the electron-like events from the 1496 days of operation of 
SuperKamiokande, published in ref. [1] and described in more detail in [46]. 
Their energy distribution is shown in fig. 2. Due to the cuts motivated by 
background rejection, the events are limited to the interval 18 — 82 MeV in 
lepton energy. 

The events could be due to neutrinos of different species. Considering only the 
dominant interaction channel for each species, we have the following possibil- 
ities (see also, e.g. [47]): 

(1) V e interacting via inverse beta decay: 

v e +p — > n + e + , (7) 

(2) u e scattering on oxygen and on electrons: 

v e + 16 -> X + e~ (8) 
v e +e~ ^>v e +e~ (9) 

(here X denotes all possible final states, X = 16 F, 16 F*, ....), 

(3) muon and tau neutrinos and antineutrinos scattering on electrons: 

v x +e~ -^v x +e~ , 

v x +e~ -*v x +e~ . (10) 

All these channels are indistinguishable in the detector because SK is not 
charge- sensitive, and so it does not discriminate between electrons and positrons. 
Moreover, since the DSNz/F is isotropic in space, one can not use angular in- 
formation to separate elastic scattering events, in contrast with an individual 
neutrino burst. 
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In the natural assumption of similar luminosities in the different species (sec. 
2.1), the process (7) largely dominates due to its larger cross section, hence 
the choice in ref. [1] to consider only this channel, obtaining a limit on the v e 
component of the DSNz/F. 

Here we study all the channels above, resulting in limits on the u e , z/ e , v x and 
v x components of the DSNz/F. We do this by considering one neutrino species 
at a time, meaning that to derive the limit on the z/ e component of the DSNz/F 
we neglected the presence of all the other species, and likewise for the other 
components of the flux. While in contrast with theory, this is an acceptable 
working assumption, as it gives the most conservative upper limits. 



> 50 n 
5 45 r 



o 40 - 




20 30 40 50 60 70 80 
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Fig. 2. The distribution in lepton energy of the SK events and of the fitted back- 
grounds. The dotted and dashed histograms below the solid one are the fitted back- 
grounds from invisible muons and atmospheric neutrinos {y e and v e ) respectively. 
The solid histogram is the sum of these two backgrounds. The dashed line above 
the solid one shows the sum of the total background and of the 90% upper limit on 
the DSNz/F signal found in [1] for the v e component of this flux. The figure is taken 
from ref. [46] , with permission. Note that the figure published in [1] is corrected for 
efficiency, i.e., it does not include the cuts that motivate the realistic efficiency in 
Eq. (12). 
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3.2 The analysis 



For each of the channels given in sec. 3.1, we performed a x 2 analysis following 
the procedure used by the SK collaboration in [1,46], to which we refer for 
details. Briefly, the steps of the analysis are as follows: 

(1) We consider the events due to the DSNi/F (signal) and those due to two 
sources of background. These are invisible muonQ and the flux of atmo- 
spheric z/ e 's and 9 e 's (figure 2). For each of these three components, we 
hold the spectral shape (taken from [46]) fixed and take the normaliza- 
tion as fit parameter. Therefore, we have three parameters, one of them 
being the DSNz/F normalization, <pQ (Eq. (6)). 

(2) We calculate the x 2 \ which is a function of the three normalizations. It 
is given by the expression: 



2 _ is (-N? ata + + gg + iQ, 

X — 2.^ I data\2 i rr2 \ ' 

1=1 \ a l > + a sys 

where (3 and 7 are normalization constants for invisible muons and for 
atmospheric electron neutrinos respectively. Ai, B\ and C\ are the num- 
ber of events from the supernova flux, invisible muons and atmospheric 
neutrinos in the I bin. jV/ fa * a is the number of events in the same bin; af ata 
and a sys are the statistical and systematic errors, from [1,46]. 

(3) We marginalize over the background normalizations and so obtain a func- 
tion x 2 (0o), which depends on the signal only. 

(4) From x 2 (0o) we find the 90% CL. limit on O - This limit is immediately 
translated into two more physically meaningful limits: one on the DSNz/F 
and another on the number of events due to it. 

(5) We repeat the procedure for different spectral shapes of the diffuse su- 
pernova flux, i.e. for different values of Eq in the interval Eq = 3.8 — 8.5 
(Sec. 2). 

To predict the histogram of events due to the DSNz/F, we have used the 
experimental efficiency as in [1,46]: 



, 0.47 if E e < 34 MeV; 
<E e ) = { e " , (12) 

0.79 if E P > 34 MeV. 



2 A muon that has too low energy to produce any Cerenkov light in water is called 
an invisible muon. Its only observable effect is the Cerenkov light from the electron 
produced by its decay. Invisible muons originate from the interaction of atmospheric 
neutrinos in the vicinity of the detector. 
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and the SuperKamiokande energy resolution as given in [48]: 



£{E e ) = 0.5 MeV^ , (13) 
with E e being the positron/electron energy. 

Published nuclear cross sections were adopted. In particular, for inverse beta 
decay we followed Vissani and Strumia [49] (fig. 3). For the energy interval 
considered here, nucleon recoil is negligible, therefore the energy of the emitted 
positron simply differs by 1.29 MeV from the neutrino energy. The total cross 
section of v e scattering on oxygen was taken from [50], with the analytical 
form in [51] (fig. 3). Following the latter reference, we have assumed that the 
difference between electron energy and neutrino energy is 15 MeV in 100% 
of the caseQ For the scattering of all species on electrons, the Standard 
Model cross section was used at the lowest order in the fine structure constant 
(see e.g., [51]). Note that the difference between the energies of the incoming 
neutrino and of the scattered electron is not a constant in this case. 

The analysis was restricted to the physical region, <ft > 0, by normalizing the 
likelihood function to 1 in this semi-plane, as it is done in [1,46]. 

We found that the data are not constraining enough to have Eq as an addi- 
tional fit parameter, instead than a fixed quantity. That would mean having 
four fit variables, two for the DSNz/F and two for the background. This num- 
ber of degrees of freedom is too large considering that the DSNz/F falls rapidly 
with energy, and therefore it is constrained mainly by the twelve events in the 
first two energy bins (fig. 2). 



4 Results 

We now present the results for each of the four channels, u e , u e , u x , v x . They 
are shown in the Tables 2-5, and in the figures 4-7. For each channel, we give 
both the 90% C.L. upper limit on the number of events, N 90 %, as well as 
the corresponding bound on the neutrino flux for certain intervals of neutrino 
energy. 

While adequate overall, the exponential spectrum we use, Eq. (6), tends to 
overestimate the flux at low energy (see [43]). Therefore our upper limits for 



3 This assumption has only a phenomenological, approximated validity. Improve- 
ments on this should come from detailed cross section calculations of neutrino scat- 
tering on oxygen, and have not been published so far. 
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Fig. 3. The cross sections for the processes (7), (8) and (9), multiplied by the relevant 
numbers of scatterers in Super Kamiokande, as functions of the neutrino energy, E. 
The cross section for ^-electron scattering is integrated over the energy of the 
outgoing electron, E e , under the condition that it exceeds the detector's threshold, 
E e > 18 MeV. Note how, in spite of the smaller cross section, ^.-electron scattering 
dominates over z^ e -oxygen interaction at E ~ 20 — 30 MeV, thanks to the larger 
number of scatterers. We used the number of protons N„ = 1.5 • 10 33 . 



thresholds around 10 — 15 MeV are slightly more conservative than those that 
could be obtained with more realistic (but model-dependent) spectral shapes. 

We do not present limits on the full fluxes (i.e., integrated over all energies) 
for two reasons. The first is that below E ~ 10 MeV the energy spectrum of 
the DSNz/F becomes more model dependent and therefore any bound could 
be given only for very specific scenarios of neutrino spectra and luminosities, 
supernova rate, etc., as done in [1]. The second reason is that the study of 
the full flux is out of the reach of current and planned experiments, due to 
background at low energ}]" 4 "] 



4 Even liquid scintillator detectors [52,53], as well as the upgraded SK with the 
addition of Gadolinium [54], would be limited to a ~ 10 MeV threshold, due to 
the ineliminable background of reactor neutrinos. Lowering the threshold could be 
feasible if a detector becomes available in a region free from sources of nuclear power 
[53]. 
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Table 2 

Upper limits for the v e channel. We give the 90% C.L. bound for the number of 
events, iVg %, and the corresponding limit on the flux of neutrinos above two differ- 
ent energy thresholds. The numbers in brackets in the second column were obtained 
using the less accurate cross section adopted in the SK analysis [1]; they are shown 
for comparison. 

4-1 Results for u e 

Our results for the v e flux are summarized in Table 2 and fig. 4. They generalize 
the study of the SK collaboration [1], where only models with Eq ~ 5.1 — 5.7 
MeV were considered. We give the flux limits above two thresholds of neutrino 
energy: E thr . = 19.3 MeV and E thr = 11.3 MeV. The first corresponds to the 
energy window used in the SK analysis in [1], while the second is the expected 
threshold for SK upgraded with the addition of Gadolinium to the water [54] 
(see sec. 5). 

The 90% C.L. bound on the number of events is in the range iVg % ~ 5 — 12, 
and the corresponding upper limit on the flux above 19.3 MeV of neutrino 
energy is ~ 1.4 — 1.9 cm~ 2 s _1 . Both the bounds on the flux and that on the 
number of events depend on Eq monotonically, getting looser for larger Eq. 
This agrees with intuition: with the decrease of Eq the neutrino spectrum and 
the spectrum of the signal of observed electrons fall more rapidly with energy. 
This implies that the DSNz/F is constrained by the uncertainty on the number 
of events in the lowest one or two energy bins (see figure 2). In contrast, with 
a less steep spectrum (larger Eq) part of the signal could be in the higher 
energy bins, resulting in a looser constraint. 

For the energy spectra used in the SK analysis [1], e.g., Eq = 5.35 MeV, 
we obtained a limit that is about 25% looser than the SK result, Eq. (1): 
$ Pe (£ > 19.3 MeV) ~ 1.49 cm^s" 1 , compared to $ Pe (£ > 19.3 MeV) ~ 
1.2 cm _2 s _1 of ref. [1]. We have checked that a 17% discrepancy is due to our 
using a more precise cross section for inverse beta decay, while the remaining 
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Fig. 4. The 90% C.L. upper limit on the v e flux (in cm 2 s 1 ) above a neutrino 
energy threshold Ethr, as a function of Ethr, for different values of Eq. 

8% difference should probably be attributed to details in the data analysis of 
the SK collaboration that we do not have access to. It should be considered 
as an error associated to the specific method of analysis used. 

The upper limit on the diffuse flux of v e above 11.3 MeV varies between 
~ 4.9 — 13.4 cm~ 2 s _1 . The allowed flux is larger for smaller Eq, i.e., for steeper 
spectrum: indeed, a spectrum that falls rapidly in energy allows a large flux 
below the SK threshold, while still giving a sufficiently small number of events 
in the SK energy window. 

Figure 4 completes the description of the v e flux bounds, by showing how 
these vary with the energy interval, which was taken of the form [£^,00]. 
This figure is of guidance for experimental projects that may have different 
energy threshold than SK. Notice how around Ethr ~ 20 MeV the flux limit 
has little dependence on Eq This because E t hr ~ 20 MeV corresponds to 

5 We did not find a value of E t h r for which the dependence on Eq is null. The curves 
in fig. 4, as well as those for other channels, figs. 5-7, do not meet in one point, even 
though the figures might give such illusion due to the use of the logarithmic scale. 
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Table 3 



The same as Tab. 2 for the v e channel. 

the interval of the SK analysis, and for this reason it is the most strongly 
constrained. Due to the exponential dependence of the DSNz/F on energy, a 
small variation of the flux - within the statistical errors of SK - at 20-30 MeV 
corresponds to a large variation of the same flux at lower or higher energy, and 
this explains the strong dependence of the bounds on E for E t hr substantially 
different from 20 MeV. 



4-2 Results for v e 

The upper limits on the v e component of the DSNz/F, and on the associated 
number of events, are given in Table 3 and fig. 5. In Table 3 we present the 
limits for the v e flux above 11.3 MeV and 19.3 MeV to allow direct comparison 
with the v e channel, where these thresholds are most relevant. We also give 
results for the same energy window of the SNO search, E = 22.9 — 36.9 MeV, 
also for the purpose of comparison. 

Similarly to the u e channel, we see that larger E corresponds to looser limit on 
the number of events, which can be as high as ~ 33. The upper bound on the 
flux of v e is 1.5 • 10 3 cm'V 1 (1.5 • 10 2 curV 1 ) above 11.3 MeV (19.3 MeV). 
It is about two orders of magnitude larger than the limit on the v e flux in the 
same energy interval, due to the smaller detection cross section. Besides the 
overall scale, the v e cross section (total of both detection processes) differs from 
that of inverse beta decay in its rising faster with energy. This implies a less 
steep fall of the spectrum of the observed leptons, and therefore an increased 
possibility to have DSNz/F events distributed in the high energy bins within 
the statistical errors. That is why the limits on the numbers of events are 
generally looser for v e . The same difference in the energy dependence of the 
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Fig. 5. The same as fig. 4 for u e . 

cross section explains why in the v e channel the correlation between lower Eq 
and larger flux allowed is more pronounced. 

In the interval E = 22.9 — 36.9 MeV we find the v e flux to be smaller than 
54 cm~ 2 s -1 . Such constraint improves on the SNO result in the same interval, 
Eq. (3), and therefore represents the best direct bound to date. 

We refer to Fig. 5 for the bounds in other energy intervals; the qualitative 
features of the figure are similar to those of fig. 4. 



4-3 Results for v x and u c 



The upper bounds on N$ Q % and on the flux for v x and 9 X , shown in Ta- 
bles 4 and 5, and figures 6 and 7, are of the same order of magnitude, with 
small differences reflecting differences in the cross sections of neutrinos and 
antineutrinos. For E > 19.3 MeV (Tables 4 and 5) the limit on the flux is 
~ (1.0 - 1.35) • 10 3 cm" V 1 for u x and ~ (1.3 - 1.8) • 10 3 cm" V 1 for V x . These 
improve by about four orders of magnitude on the previous best limits from 
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£//MeV > 11. 
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1.32 • 10 4 
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3 8 

O.O 


A M 
4.01 
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1.31 • 10 3 
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o.zy • 10 
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5.32 


4.99 • 10 3 
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6.5 


5.87 


3.78 • 10 3 


1.11 • 10 3 


7.5 


6.45 


3.08 • 10 3 


1.06 • 10 3 


8.5 


7.06 


2.66 • 10 3 


1.02 • 10 3 



Table 4 

The same as Tab. 2 for u x . 



Eo/MeV 


^90% 


limit on (cm 2 s 1 ) 






E/MeV > 11.3 


E/MeV > 19.3 


3.5 


4.38 


1.74 • 10 4 


1.77- 10 3 


3.8 


4.51 


1.41 • 10 4 


1.72 • 10 3 


5.35 


5.22 


6.81 • 10 3 


1.53 • 10 3 


5.5 


5.30 


6.48 • 10 3 


1.51 • 10 3 


6.5 


5.82 


4.88 • 10 3 


1.43 • 10 3 


7.5 


6.40 


3.96 • 10 3 


1.36 • 10 3 


8.5 


7.00 


3.35 • 10 3 


1.31 • 10 3 



Table 5 

The same as Tab. 2 for v. 



LSD, Eq. (4), which refer to the interval 20 - 100 MeV. 

Due to the rough equipartition of energy among the six neutrino species pre- 
dicted by theory, it is reasonable to expect the u x and u x components of the 
DSNz/F to be comparable to the electron flavor ones, that are more strongly 
constrained. Therefore, the direct limits we found here for v x and u x , while 
being the best available, are probably far from realistic values of these fluxes. 
They are nevertheless important for their observational, model-independent 
character, and as a reference for future, more sensitive, searches. 
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Fig. 6. The same as fig. 4 for v x . 
4-4 Interpretation of the results: luminosity bounds 



What constraints can be drawn from our results on the neutrino emission from 
a supernova? A full address of this question requires a detailed analysis of all 
the uncertainties, which we do not perform here. However, as an example we 
study what can be concluded on the neutrino luminosity if the supernova rate 
is known, which is likely to be the case in the near future, after data from 
SNAP [55] and JWST [56] become available. 

To extract the luminosity in a given flavor of an individual neutrino burst, 
one needs to go beyond the simple parametrization of the diffuse flux in Eq. 
(6). We use the analytical description by Lunardini [43], which is valid within 
a few per cent above the SK threshold and is obtained from Eq. (5) with a 
power law form of Rsn and a neutrino spectrum at the exit of the star (after 
oscillations) of the form of an exponential times a power law [21]. It reads: 



Rsn(0)Lc (E_ 
H T(2 + a)Ei Uo 



a—ri—l 



»7 + l, 



E 
En 



T} + 1,(1 + Zr, 



I E 

En 



(14) 
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Fig. 7. The same as fig. 4 for v x . 

where L is the total energy (in the individual flavor) of the neutrino burst 
leaving the star, and a ~ 2 — 5 is a parameter describing its spectral shape 
[43]. The average energy of the neutrinos leaving the star is (E) ~ E (l + a). 
We have r\ = a+j3—3Q m /2, where (5 describes the cosmo logical supernova rate: 
R SN ~ (1 + z)P . z max ~ 1 is a break in the redshift distribution of supernovae. 
We fix R SN to be the best piecewise fit to the data found by Beacom and 
Hopkins [40]: R SN ^ 2.0 • 10~ 4 yr^Mpc" 3 , (5 = 3.28 and z max = 1.0. While 
Eqs. (6) and (14) differ only minimally in the energy dependence [43], the 
second is more adequate for our purpose because in it the dependence on the 
different parameters is explicit. 



For every value of Eq and a, the combination of Eq. (14) and of our flux 
limits (Tables 2-5) gives an upper limit on the integrated luminosity L, which 
can then be compared to the theoretical expectation of L - 1 1 ) J 
neutrino flavor. 



5 • 10 ergs per 



Results are shown in fig. 8, where we have imposed 9 < (E)/MeV < 24 as 
naturalness constraint (solid curves). We see that, for E <: 5.5 MeV, the 
upper limit on the V e luminosity approaches, or even crosses, the theoretically 
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interesting range, an intriguing fact that has been already observed in [57]. 
Softer neutrino spectra, however, allow neutrino the neutrino luminosity to be 
several times larger than the theoretical central value. This is expected, since 
softer neutrino spectrum means smaller fraction of the total energy above the 
SK threshold. A further loosening of the luminosity constraint is expected if 
the uncertainty on the supernova rate is included [15]. 

The luminosity bounds for the other neutrino species are far from the theoret- 
ical range, reflecting the weaker flux limits. They could be useful to constrain 
the presence hidden sources that could add to the expected supernova neutrino 
flux. 



5 Summary and discussion 

We have analyzed the data of the search for the diffuse supernova neutrino 
flux published by the SK collaboration in 2003 [1]. Our work extends the 
analysis of collaboration itself in a number of ways. First, it considers all the 
relevant channels, one at a time (i.e., with all the signal attributed to a certain 
channel), resulting in bounds on the fluxes of all the neutrino species. Second, 
it takes into account a wider range of neutrino spectra, including the softer 
ones that have been suggested in the literature recently (Table 1). It also uses 
a more precise cross section for inverse beta decay, that rises more slowly with 
energy and therefore implies a looser bound on the v e component of the flux. 
Results have been given for several different energy intervals motivated by 
near future experiments. 

A first result is that the flux of v e is bound, at 90% CL, to be smaller than 
1.4-1.9 cm-V 1 for E > 19.3 MeV and than 4.9-13.4 cm^s' 1 above a 11.3 
MeV threshold. The flux intervals account for the variation with the neutrino 
energy spectrum, and the discrepancy with the SK limit, Eq. (1), is explained 
by the different cross section used, up to a minor 8% difference. 

The limits in the v e channel improve on the previous limit from SNO, Eq. 
(3): at 90% C.L. the v e flux must be smaller than 54 cm~ 2 s _1 in the interval 
22.9 — 36.9 MeV even for the softest neutrino spectrum, which represents the 
most optimistic case. The flux is allowed to be as large as (9(10 2 ) cm~ 2 s~ 1 
(C(10 3 ) cm-V 1 ) above 19.3 (11.3) MeV energy. 

The upper limits on the non-electron species are in the range of O(10 3 ) cm~ 2 s _1 
to (9(10 4 ) cm _2 s _1 for the same two thresholds, and improve by several orders 
of magnitude over the previous bounds from LSD, Eq. (4). 

Our results represents the most one can get from this generation of detectors 
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1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 5.0 

a a 



Fig. 8. 90% C.L. upper limits on the neutrino luminosity (in units of Lq = 10 52 ergs) 
in the different species at the exit of the star, as a function of the "shape" parameter 
a. The curves, from upper to lower, correspond to Eo/MeV = 3.5,5.5,6.5,7.5,8.5 
in every panel. The solid curves were obtained by imposing a naturalness con- 
dition on the average energy of the neutrino spectrum at the exit of the star: 
9 < (E) /MeV < 24. Notice the large difference in the vertical scales between the 
different panels. 
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Table 6 



Flux limits from SK at 90% C.L. for the energy intervals and channels that are 
relevant to the four most discussed future DSNz/F detectors (see text for details). 
We distinguish between generic designs and detailed studies that include estimates 
of the background at specific sites (see Table 8 in [62]). The intervals in the flux 
correspond to varying the neutrino spectrum in the range shown in Table 1. 

on the diffuse flux, up to minor improvements that might come from increased 
statistics at SuperKamiokande (see the preliminary results in [58]). Consider- 
ing that theoretical predictions of the DSNz/F in each neutrino species range 
from touching the current SK limit for v e down to values 20 times smaller (Tab. 
1), it is clear that SK might not be able to see the DSNz/F in the v e channel 
and almost certainly will never detect it in other channels. This strongly rein- 
forces the already strong case for new detectors with enhanced sensitivity to 
the DSNz/F. 

With this in mind, we have given constraints on the flux in different energy in- 
tervals that could be relevant for future searches. They represent the minimum 
sensitivity that these searches should have to improve on the current status. 
To summarize, and also to add detail, in Table 6 we give the SK bound on the 
DSNz/F for the specific setup (neutrino species and energy interval) of the four 
different concepts that are most discussed for the future. These are a Mega- 
ton water Cerenkov detector [59,60,61], a 50 kt liquid scintillator experiment 
[52,53], a liquid Argon chamber [64,65] and the SK tank with Gadolinium 
trichloride dissolved in it [54]. For some projects, we consider both generic 
and site-specific designs, for which the background, and therefore the energy 
window of sensitivity, have been calculated in detaif^J From the Table one 

6 The energy intervals for the site-specific designs are nevertheless tentative, be- 
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can immediately see how crucial it is to have a detection threshold as low as 
possible, since the DSNVF falls exponentially with energy. 

The SK bounds in Table 6 will certainly be a practical benchmark for the 
planning of future detectors, as they allow to give a quantitative estimate of 
how much a given design will improve over SK. Moreover, our bounds will be 
an important component of a global analysis of multiple data sets once new 
data are available from different experiments. 

We conclude with the comment that the detection of the diffuse flux of su- 
pernova neutrinos is an opportunity that the underground science community 
can not afford to miss, because of its implications for fundamental physics 
and astrophysics/cosmology. Since the flux is continuous in time, the barrier 
to its observation is only technological. Once the barrier is overcome, the ex- 
perimental study of supernova neutrinos will transition from the realm of rare 
event to that of regular, ongoing data taking and analysis. Our results repre- 
sent the status of the art of current searches and the ultimate sensitivity - up 
to minor improvements - of the existing detectors. It marks where the present 
generation of experiments leaves off and the new one will take on, in the next 
decade or so, to bring this field to maturity. 
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